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Abstract 

Simulation  of  semiconductor  diode  laser  performance  involves  interaction  be¬ 
tween  multiple  physics  domains.  This  report  presents  the  governing  equa¬ 
tions  and  finite-element  (FE)  implementation  of  such  a  simulation,  includ¬ 
ing  coupled  calculation  of  electronic  band  structure,  determined  by  solving 
Schrodinger’s  equation;  carrier  distribution  and  potential  in  the  active  region, 
solved  using  a  drift-diffusion  model;  electrostatic  potential  in  the  ohmic  re¬ 
gions,  governed  by  Laplace’s  equation;  temperature  distribution,  governed 
by  Fourier’s  law  and  the  heat  equation;  and  light  propagation,  described  by 
Maxwell’s  equations.  Preliminary  simulation  results  are  provided  for  stripe 
lasers  based  on  the  GaAs/AlGaAs  material  system,  including  asymmetri¬ 
cally  waveguided  devices.  Important  implementation  issues  are  discussed, 
particularly  computation  across  multiple  length  scales  and  formulation  of 
appropriate  boundary  conditions.  Simulation  results  are  compared  to  pre¬ 
vious  calculations  reported  in  the  literature.  Advancement  of  the  technical 
component  is  paralleled  by  development  of  a  graphical  user  interface  for  ease 
of  use  by  non-experts  in  FE  methods. 


1  Introduction:  Multiphysics  Modeling 

At  the  most  fundamental  level,  the  generation  and  characteristics  of  the 
light  produced  by  a  semiconductor  laser  is  due  to  the  quantum  mechanical 
interaction  of  electrons,  holes,  photons,  and  phonons.  Therefore,  electron¬ 
ics,  optics,  and  temperature  ultimately  determine  the  performance  of  a  laser 
diode.  The  magnitude  of  injection  current  needed  to  operate  a  high  powered 
laser  diode  results  in  high  temperatures,  reducing  the  external  efficiency  and 
operating  lifetime.  A  comprehensive  simulation  models  the  coupling  of  op¬ 
tics,  temperature,  and  electronics  in  a  self-consistent  manner.  However,  such 
comprehensive  models  often  face  numerous  computational  problems.  One 
such  problem  is  due  to  the  multiple  length  scales  in  a  laser  structure.  The 
thicknesses  of  the  layers  in  a  given  diode  vary  between  several  orders  of  mag¬ 
nitude.  For  example,  the  thickness  of  the  substrate  is  typically  on  the  order  of 
a  hundred  microns,  whereas  the  thickness  of  the  quantum  well  is  less  than  ten 
nanometers.  In  thick  layers,  classical  characteristics  dominate  the  behavior 
of  carriers.  However,  in  the  quantum  well,  the  band  structure  becomes  quan¬ 
tized  in  energy.  Further,  the  radical  differences  in  the  layer  thicknesses  leads 
to  computationally  taxing  simulations,  owing  to  the  mesh  densities  required 
for  an  accurate  solution.  Therefore,  a  useful  simulation  must  capture  the 
essential  physical  phenomena  and  simultaneously  be  numerically  tractable 
for  the  computational  hardware  available.  We  have  developed  a  preliminary 
graphical  user  interface  that  makes  use  of  the  Matlab  scripting  options  in 
COMSOL  Multiphysics.  This  allows  the  user  to  perform  simulations  with¬ 
out  expertise  in  the  finite-element  method  and  still  maintain  control  over  the 
device  geometry  and  composition. 

2  Electrothermal  Model 

The  electrothermal  model  solves  for  the  potential  and  temperature  fields, 
models  the  electron  and  hole  transport.  Ultimately,  the  current- voltage  (IV) 
characteristics  of  the  diode  are  given  by  the  transport  model.  Because  of 
the  interaction  between  the  drift  and  diffusion  of  carriers,  the  IV  character¬ 
istics  are  non-Ohmic  or  nonlinear.  Joule,  or  resistive,  heating  depends  on 
both  the  current  density  and  the  electric  held;  and  the  Fermi-Dirac  distribu¬ 
tions  that  determine  the  nonequilibrium  carrier  densities  depend  on  tempera¬ 
ture.  Therefore,  the  temperature  held  is  strongly  coupled  to  the  electron-hole 
transport  model.  The  default  simulation  structure  throughout  this  report  is 
similar  to  the  one  presented  in  [17J.  A  contour  plot  of  the  temperature  over 
a  mirrored  facet  of  the  laser  is  shown  in  Fig.  2, 
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Figure  1:  Schematic  illustrating  the  coupling  between  physical  domains  in 
semiconductor  lasers. 
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Figure  2:  Contour  plot  of  temperature  over  a  mirrored  facet  of  a  high  powered 
diode  laser. 


Optical  heating  is  due  to  the  absorption  of  photons  by  free  carriers,  and 
therefore  depends  on  the  photon  and  free  carrier  densities  |21j.  Optical  heat¬ 
ing  is  not  modeled  in  the  present  work,  but  is  particularly  important  in  long 
wavelength  lasers  emitting  around  1.55  /im  [2J. 

2.1  Drift-Diffusion  Model 

The  drift-diffusion  equations  are  used  to  model  the  transport  of  electrons 
and  holes  through  the  device.  The  nomenclature  of  the  equations  results 
from  the  phenomena  they  model:  carrier  drift,  due  to  the  electric  field;  and 
carrier  diffusion,  due  to  concentration  gradients.  The  electrostatic  potential 
is  determined  by  the  solution  to  Poisson’s  equation, 

V  •  (—eVip)  =  q  (p  —  n  +  N)  (1) 

where  ip  is  the  electrostatic  potential,  e  =  es£0  is  the  permittivity,  q  the 
electron  charge,  p  the  hole  concentration,  n  the  electron  concentration,  and 
IV  =  Np  —  Nj  the  doping  profile.  The  model  assumes  the  full  ionization  of 
dopants. 

Further,  the  drift-diffusion  equations  solve  for  the  electron  and  hole  con¬ 
tributions  to  the  current  density.  The  current  densities  are  determined  by 
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continuity  equations  that  balance  the  flux  of  carriers  through  a  control  vol¬ 
ume  with  the  rate  at  which  carriers  recombine.  For  electrons,  that  is, 


dn 

dt 


q  ^ 


(2) 


where  jn  is  the  electron  contribution  to  the  current  density,  and  R  is  the 
sum  of  the  net  recombination  rates,  which  are  defined  as  the  recombination 
rates  minus  the  generation  rates.  An  analogous  equation  holds  for  the  hole 
flux, 

(3) 

where  jp  is  the  hole  contribution  to  the  hole  current  density.  The  total  current 
density  is  given  by  the  j  =  jn  + j  .  The  electron  and  hole  contributions  to  the 
current  density  are  given  by  the  quasi-Fermi  formulation  of  the  drift- diffusion 
equations  |29j,  [14J.  The  electron  current  density  is  calculated  by, 


jn  =  -qnUnVtfn- 


(4) 


Similarly,  the  hole  current  density  is, 

j  P  =  -gp/ipVyy  (5) 

The  continuity  equations  are  implemented  in  steady  state,  so  that  the  time 
derivatives  are  zero.  The  quasi-Fermi  potentials  are  related  to  the  quasi- 
Fermi  energies  as  Epc  =  —  qipn  and  Epv  =  —q^pP-  The  electron  and  hole 
mobilities  are  given  by  pn  and  pp,  respectively.  The  calculation  of  carrier 
mobility  varies  in  complexity,  ranging  from  simple,  constant,  low-field  mobil¬ 
ities;  to  computationally-taxing  molecular  dynamics  simulations.  Since  the 
current  densities  are  directly  proportional  to  the  mobility,  is  important  that 
it  be  accurately  determined.  Additionally,  the  mobility  must  be  computa¬ 
tionally  tractable  in  a  reasonable  time,  and  should  reflect  the  dependence 
of  mobility  on  temperature.  The  mobilities  in  the  present  work  are  calcu¬ 
lated  using  the  semi-empirical,  Caughey-Thomas-like  parameters  given  in 
|27j.  The  mobility  model  gives  the  dependence  of  temperature  and  doping, 
and  is  well-suited  for  device  simulation. 


2.1.1  Carrier  Concentrations 


In  general,  the  carrier  concentrations  can  be  determined  by  integrating  the 
density  of  states  and  occupation  functions  over  all  energies  in  a  energy  band, 
such  as  _ 


n  =  /  pc  (E)  fc  (E)  dE 
J  Ec 


(6) 
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for  the  conduction  band,  where  Ec  is  the  conduction  band  minimum,  pc  is 
the  density  of  states  function,  and  fc  is  the  Fermi-Dirac  function  for  the 
conduction  band.  A  significant  simplification  can  be  made  by  assuming  non¬ 
degenerate  doping  levels  |5j.  This  assumption  is  justified  if  the  doping  is 
such  that  the  equilibrium  Fermi  energy  EF  is  at  least  3 kBT  less  than  the 
conduction  band  minimum,  and  at  least  3 kBT  greater  than  the  valence  band 
maximum.  Under  these  conditions,  the  Fermi-Dirac  distribution  approaches 
the  Maxwell- Boltzmann  distribution,  and  the  bulk  carrier  density  can  be 
calculated  from  the  quasi-Fermi  energies  without  the  need  for  integration,  as 

n  =  rii  exp  [q  -  (pn )  /  (kBT)]  =  Nc  exp  [(EFc  -  Ec )  /  (kBT)\ ,  (7) 

and 

P  =  ni  exp  [q  {(pp  -  ip)  /  (. kBT )]  =  Nv  exp  [(Ev  -  EFv)  /  (kBT)\ ,  (8) 

where  Nc  and  Nv  are  the  effective  density  of  states  for  the  conduction  and 
valence  bands,  respectively.  These  are  given  by  Nc  =  2  f 2irm*chBT/h 2]3^2  and 

Nv  —  2  [2nmlkBT / h2]3^2 ,  where  m*v  =  ( m *^j2  +  m ^  for  light  and  heavy 
holes  [14J,  [5j.  The  equalities  above,  wherein  the  carrier  density  is  a  function 
of  the  quasi-Fermi  energy,  are  useful  in  the  calculation  of  the  recombination 
rates,  which  depend  on  the  quasi-Fermi  energies  relative  to  the  band  extrema. 

The  above  carrier  density  calculation  is  valid  in  the  bulk.  However,  quan¬ 
tization  of  the  energy  bands  leads  to  a  different  density  of  states  function  for 
carriers  with  energy  less  than  the  quantum  well  barrier  energy.  Additionally, 
some  carriers  may  have  energy  greater  than  the  quantum  well  barrier;  these 
carriers  obey  the  bulk  carrier  density  equations  (7)  and  (8).  The  total  carrier 
density  calculated  at  the  quantum  well  is  the  sum  of  these  two.  The  density 
of  bound  carriers  in  the  quantum  well  is  given  by  [18J ,  [7j  as 

Abound  =  In  I1  +  exp  [-  ( Ecn  -  EFc )  /  (kBT)\ }  ,  (9) 

7T  Ti  da  n 

where  the  summation  is  over  the  bound  states  computed  in  Section  5,  An 
an  analogous  equation  is  used  for  bound  holes. 

2.2  Conduction  Heat  Transfer  Model 

The  temperature  held  is  computed  using  Fourier’s  law  and  the  heat  equation. 
The  present  model  considers  only  Joule  heating,  although  optical  heating  has 
proven  to  be  significant  in  some  material  systems,  particularly  InGaAsP- 
based,  long  wavelength  lasers  [21J.  The  heat  equation  is, 
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V  ■  (-kvt)  =  (10) 

where  k  is  the  thermal  conductivity,  T  is  the  lattice  temperature.  The  Joule 
heating  term  is  calculated  using, 

Qj'oule  =  -J  •  (11) 


2.3  Boundary  Conditions  for  Electrothermal  Model 

The  boundary  conditions  for  the  electrothermal  model  are  shown  in  Fig. 
3.  Far  from  the  depletion  regions  at  the  device  contacts,  it  is  reasonable 
to  assume  that  the  carrier  concentrations  relax  to  their  equilibrium  values 
([4],  pg  313;  [lj).  The  potential  will  also  be  at  its  equilibrium  value  on 
the  n-contact,  calculated  by  eq.  (12).  If  the  device  is  in  forward  bias  with 
applied  voltage  Va,  the  p-contact  will  be  at  a  potential  of  if>0  +  Va.  In  the 
framework  of  quasi-Fermi  levels,  this  implies  that  the  quasi-Fermi  energies  far 
from  the  junction  return  to  the  equilibrium  Fermi  level  Ep,  which  describes 
the  equilibrium  carrier  statistics.  If  the  Fermi  potential  (tpp  =  — Ep/q )  is 
used  as  the  reference,  it  may  be  shown  that  the  n-contact  uses  the  Dirichlct 
conditions  <pn  =  0 ,  (pp  =  0,  and  if}  =  ^o-  The  p-doped  contact  uses  the 
Dirichlet  conditions  ipn  =  Va,  (pp  =  Va,  and  if}  —  ij}, o  +  Va.  The  metal- 
semiconductor  interface  physics  are  not  modeled  because  the  semiconductor 
materials  near  the  n-  and  p-contacts  are  heavily  doped  (n+  and  p+)  to  prevent 
a  large  voltage  drop  across  this  barrier  (pg.  175  of  |9J).  It  is  also  worth  noting 
that  the  applied  voltage  must  be  gradually  increased  from  zero  in  increments 
of  approximately  ksT/q  m  0.026  V  (at  room  temperature)  until  the  desired 
potential  is  reached  [29J. 

The  equilibrium  potential  distribution  is  a  result  of  the  diffusion  of  carri¬ 
ers  across  the  p-n  junction,  which  sets  up  an  electric  held.  The  equilibrium 
potential  distribution  is  computed  with  [lj , 


V>o=  —  (In 
Q 


nm 

rij 


+  In 


Ppo 

rij 


with  equilibrium  electron  concentration  in  the  n-doped  region, 

|  IV  |  +  yj. N 2  +  4  nf 


nm  — 


and  equilibrium  hole  concentration  in  the  p-doped  region, 

|  IV  |  +  yj N 2  +  4  nf 


Ppo  = 


(12) 


(13) 


(14) 
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n-contact  (thermally  insulated,  electrical  ground) 


=  0  <p„  =  0 

n  ■  (kV T)  =  0  <pp  =  0 


n  -  V4>  =  0 

n  •jp  =  ° 

«  •jn  =  0 

n  ■  (k 7T)  = 
(thermally, 
electrically 
insulated) 


n  •  =  0 

n  -jp  =  0 

n  '  J  n  =  0 
n  ■  (kVT)  =  0 
(thermally, 
electrically 
insulated) 


Figure  3:  Boundary  conditions  for  the  electothermal  model.  The  p-contact  is 
attached  to  an  ideal  heat  sink  at  T  =  Ths,  and  at  drive  voltage  ip  —  ip0  +  Va. 
The  n-contact  is  thermally  insulated,  and  serves  as  electrical  ground. 


The  thermal  boundary  condition  over  the  top  contact  has  been  investi¬ 
gated  experimentally  |28j  and  it  was  determined  that  due  to  the  presence 
of  a  metal  electrical  contact,  the  boundary  should  be  modeled  as  isother¬ 
mal,  rather  than  convective  or  insulated.  This  is  naturally  accounted  for  in 
our  model,  since  the  metal  contact  included  in  the  simulation  has  high  ther¬ 
mal  conductivity  and  therefore  facilitates  a  uniform  heat  distribution.  All 
other  boundaries  are  thermally  insulated,  with  the  exception  of  the  boundary 
in  contact  with  the  heat  sink,  which  uses  the  Dirichlet  boundary  condition 
T  =  Ths  to  account  for  the  heat  sink  temperature.  This  models  the  heat  sink 
as  ideal. 


3  Recombination  Mechanisms 

The  process  by  which  an  electron  in  the  conduction  band  falls  into  an  empty 
state  in  the  valence  band  is  known  as  recombination.  The  potential  energy 
lost  by  the  electron  during  recombination  is  manifested  in  another  form.  In 
the  event  that  the  energy  is  transformed  to  light,  by  the  emission  of  a  photon, 
the  process  is  known  as  radiative  recombination.  Stimulated  and  spontaneous 
recombination  are  radiative  recombination  mechanisms.  All  other  processes 
are  loosely  termed  nonradiative  recombination.  Nonradiative  recombination 
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can  result  in  the  creation  of  phonons,  through  Shockley-Read-Hall  recombi¬ 
nation;  or  the  energy  can  be  kinetic  energy  imparted  to  an  electron  or  hole, 
by  way  of  Auger  recombination  (pg.  193-194  of  |5j). 

3.1  Stimulated  Recombination 

The  key  feature  in  the  operation  of  a  semiconductor  is  stimulated  emission  via 
optical  gain.  Optical  gain  is  achieved  by  the  condition  of  population  inver¬ 
sion,  wherein  electrons  are  stimulated  to  the  conduction  band,  leaving  holes 
in  the  valence  band.  Population  inversion  occurs  when  the  potential  across 
the  quantum  well  leads  to  a  quasi-Fermi  level  separation  greater  than  the 
band  gap  |9j,  |7J.  When  population  inversion  has  been  achieved,  optical  gain 
occurs  by  the  interaction  of  the  electromagnetic  field  with  the  electron-hole 
pairs,  leading  to  recombination  which  emits  a  photon  of  the  same  frequency 
as  the  perturbing  electromagnetic  field.  Optical  gain  gives  the  increase  per 
unit  length  in  the  number  of  photons,  and  thus  its  calculation  depends  on 
the  quasi-Fermi  energy  separation  determined  by  the  drift-diffusion  model. 
The  stimulated  emission  rate  (in  units  of  per  unit  volume  per  unit  time)  is 
related  to  the  optical  gain  by, 


Rst  VgfjNph  (15) 

where  vg  is  the  group  velocity  and  Nph  the  photon  density.  The  calculation 
of  these  parameters  is  discussed  in  the  subsequent  subsections. 


3.1.1  Optical  Gain 

Optical  gain  is  defined  as  the  proportional  increase  in  the  photon  density  per 
unit  length  of  active  material.  Parabolic  band  shapes  and  ^-selection  rules 
are  assumed.  The  expression  for  the  optical  gain  as  a  function  of  the  photon 
energy  is  [3j,  [30],  ch.  1  and  2  of  [18J, 


g{Ei)  = 


nq2hn, 


g,eS 


£0mlc0rfEi 


EE 


Eq-\-ECn-\-Ej ■ 


\Mqw\2  Pred,j  (E)  ( fcn,j  ~  fvn,j )  £  (E)  dE . 


(16) 

Here,  Co  is  the  speed  of  light  in  a  vacuum,  m,Q  is  the  free  electron  mass,  nr 
is  the  refractive  index,  ng>e g  is  the  group  effective  refractive  index,  Et  is  the 
photon  energy,  and  j  =  hh,  Ih  for  heavy  holes  and  light  holes,  respectively. 
|Mqh/|2  is  the  squared  transition  matrix  element  for  a  quantum  well  given 


by  |30j, 


\Mqw\2  =  ^  \Mavg\2 


1  + 


Ecn  + 


red,j 


E 

J-'cn 

(E  -  Eg-  Ecn  -  E 


jnj 


(17) 


with 

i  t\ ,r  ,2  _  m0(m0  1  \  Eg  (Eg  +  A) 

1 avgl  ~  6  W*  )  Eg  +  2A/3  ' 


(18) 


m*  is  a  modified  effective  mass  that  takes  account  for  higher  and  lower  con¬ 
duction  bands  (taken  to  be  0.053mo  in  GaAs),  A  is  the  spin-orbit  splitting 
energy,  and  Eg  is  the  band  gap.  The  derivation  of  \Mavg\2  assumes  that  a 
single  factor  of  two  is  included  in  the  reduced  density  of  states  to  account  for 
spin  degeneracy  [30j.  The  reduced  density  of  states  function,  including  the 
factor  of  two  for  spin  degeneracy,  is  given  by  [18J,  [20J, 


Pred,j  (A^) 


^‘Kedj 

h2da 


Ti  (E  —  Eg  —  Ecn 


(19) 


where  7i  is  the  Heaviside  unit  step  function.  The  reduced  effective  mass  is 
given  by, 

c  J  (20) 


771  rpr]  a  —  j 

J  rn*  +  m* 

and  the  modified  Fermi  distribution  for  conduction  band  states  is  given  by, 

-l 


fcnj  =  { 1  +  exp 


Ecn  + 


m: 


red,j 


mi 


{E  —  Eg  —  Ecn  —  Ejn)  —  Epc  /  {ksT) 


(21) 


and  the  modified  Fermi  distribution  for  valence  band  states  is  given  by, 


fvnj  =  { 1  +  exp 


~E„ 


m: 


redd 


rn:; 


(E  —  Eg  —  Ecn  —  Ejn )  —  EFv  ]  /  (ksT) 


-i 


(22) 

Lastly,  the  spectral  broadening  function  is  a  simple  Lorentzian  distribution, 
given  by  [18J  as, 

C  (E)  =  - - 2,  (23) 

1  ;  vr  {E-Etf  +  ih/Tn)2'  K  ’ 


where  Ttn  is  the  intraband  relaxation  time,  typically  taken  to  be  0.1  ps  |3J, 
[18J.  Intraband  relaxation  or  spectral  broadening  takes  account  for  the  un¬ 
certainty  in  the  energy  states  [7j.  The  material  gain  spectrum  is  plotted  as 
as  function  of  photon  energy  in  Fig.  4, 
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Figure  4:  Material  gain  as  a  function  of  photon  energy  for  8  nm  GaAs  single 
quantum  well  device  with  a  carrier  concentration  of  1019  cm-3. 


3.1.2  Photon  Group  Velocity 

The  photon  group  velocity  is  defined  as  |9J,  [18J 

c0 

vg  = - , 

ng,eS 


(24) 


where  Co  is  the  speed  of  light  in  a  vacuum,  and  ngjeg  is  the  group  effective 
refractive  index,  which  is  given  approximately  by, 


chh,e ff  . 
<7A0 


(25) 


The  group  refractive  index  is  found  by  using  a  forward  difference  approxi¬ 
mation  to  dnr  ,eff/  <9A0. 


3.1.3  Photon  Density 

The  number  of  photons  per  unit  volume  in  the  cavity  is  determined  by  an 
energy  balance  written  for  the  optical  cavity.  The  equation  balances  the 
contributions  of  stimulated  emission,  spontaneous  emission,  and  the  rate  of 
photon  loss  via  the  mirror  and  internal  losses.  The  steady-state  balance 
equation  for  the  photon  density  in  the  quantum  well  is  |9j,  [7],  [14], 
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VgfjNph  +  R'sp  ~  Y  (ai  +  ‘lL  ln  7^7^  )  Nph  =  °-  (26) 

The  term  R'sp  gives  the  spontaneous  emission  emitted  into  the  lasing 
mode,  T  is  the  optical  confinement  factor,  cq  is  the  internal  losses  (taken  to 
be  cq  =  5  cm-1  [TJ ) ,  and  R±,  R2  are  the  mirror  refectivities.  Solving  equation 
26  for  Nph  yields, 


where  the  spontaneous  emission  into  the  lasing  mode  is  given  by  |7j, 

r-,/  _  Fvggnsp 

n sp  ~  y 

where  V  is  the  quantum  well  volume,  with  the  population  inversion  factor 
given  by, 

nsp  =  {1  -  exp  [(Ei  -  A EF)  /  ( kBT)}}~ 1 .  (29) 

where  A Ep  =  EFc  —  EFv  is  the  quasi- Fermi  level  separation. 

3.2  Spontaneous  Recombination 

For  injection  currents  below  the  threshold  current,  Rsp  3>  Rst ■  The  emitted 
light  is  of  many  polarizations  and  frequencies.  Within  the  electrothermal 
drift-diffusion  simulation,  spontaneous  emission  is  not  modeled.  However,  the 
optical  spectrum  of  spontaneous  emission  can  be  determined  by  a  calculation 
similar  to  that  of  the  optical  gain  spectrum.  The  expression  used  to  calculate 
the  spontaneous  emission  rate  is  [18J , 


(28) 


r,r  (Ei)  =  -  Y.  Y.  I  AW 2  Predj  ( E )  fcnJ  (1  -  fvnJ) ,  (30) 

TTUIqU  C0e 0  n  j 

where  nr  is  the  refractive  index  of  the  quantum  well  material,  Et  is  the  photon 
energy,  and  \Mavg\2  is  the  magnitude  of  the  energy-independent  momentum 
matrix  element,  given  by  equation  (18). 

The  total  spontaneous  emission  rate  is  found  by  integrating  over  all  pho¬ 
ton  energies  |6j  so  that, 

COO 

Rsp  =  /  t Sp  (E)  dE,  (31) 

Jo 
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Figure  5:  Spontaneous  emission  rate  as  a  function  of  photon  energy  for  8  nm 
GaAs  single  quantum  well  device  with  a  carrier  concentration  of  1019  cnr3. 


ft  is  not  necessary  to  convolve  the  spontaneous  emission  function  with  a 
broadening  function  to  determine  the  total  rate,  but  it  is  appropriate  in  a  cal¬ 
culation  of  the  spontaneous  emission  spectrum  [18J,  |30j,  |6j.  The  calculated 
spontaneous  emission  spectrum  is  shown  in  Fig.  5, 


3.3  Nonradiative  Recombination 


Two  types  of  nonradiative  recombination  present  in  a  laser  diode  are  Shockley- 
Read-Hall  (SRH)  recombination  and  Auger  recombination.  SRH  recombina¬ 
tion  is  due  to  the  presence  of  point  defects,  and  the  rate  at  which  carriers 
recombine  due  to  SRH  recombination  is  given  by,  [lj , 


Rsrh  — 


np  —  n | 

Tp  (n  +  Hi)  +Tn(p  +  Hi) 


(32) 


4  Optical  Model 

The  optical  model  calculates  the  optical  mode  shape.  The  optical  mode  shape 
determines  the  optical  confinement  factor  T,  which  describes  the  fraction  of 
light  that  is  confined  to  the  quantum  well.  The  remaining  fraction  is  absorbed 
by  the  waveguiding  layers  and  leads  to  optical  heating.  The  optical  model 
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reduces  the  diode  geometry  to  the  optical  region,  which  we  define  as  the 
portion  of  the  waveguiding  structure  where  a  non-negligible  portion  of  light 
is  confined.  This  includes  the  n-  and  p-doped  cores  and  cladding  layers,  and 
the  quantum  well. 


4.1  Calculation  of  Optical  Mode 

The  optical  mode  shape  is  determined  by  the  solution  of  the  Helmholtz  equa¬ 
tion,  which  is  derived  from  Maxwell’s  equations, 


V2E0  +  &Q  (nr  —  nr,eS )  Eo  —  0,  (33) 

where  E0  =  E0  (x,  y )  is  the  optical  mode  shape,  nr^  is  an  eigenvalue,  ko  = 
27t/Ao  is  the  wave  number  in  a  vacuum,  and  nr  is  the  complex  refractive 
index,  which  contains  information  about  the  gain  or  absorption  coefficient, 

n?  =  +  lw)  (34) 

in  the  quantum  well,  where  fir  is  the  real  part  of  the  refractive  index,  calcu¬ 
lated  by  the  model  presented  in  [10J  for  AlzGai_xAs.  In  the  passive  layers, 
the  substitution  g  — »  —a  is  made  to  account  for  absorption.  This  link  is 
not  yet  implemented  in  the  present  work,  because  the  temperature  must  be 
mapped  from  the  electrothermal  simulation  to  the  optical  model  to  give  the 
temperature  dependence  of  nr. 

Since  ko  depends  on  Ao,  Schrodinger’s  equation  must  first  be  solved  to 
determine  the  wavelength.  The  Helmholtz  equation  is  solved  in  the  x-y  plane 
of  the  laser,  with  the  wave  propagating  in  the  ^-direction.  The  boundary 
conditions  used  to  solve  the  optical  simulation  assume  that  the  mode  shape 
goes  to  zero  at  the  boundaries  of  the  optical  region.  The  simulation  uses  a  y 
cross-section  of  the  fundamental  waveguide  mode  profile  and  normalizes  the 
profile  such  that, 


/OO 

Eo  (y)dy  =  1  (35) 

-OO 

A  normalized  optical  mode  distribution  for  a  symmetric  waveguide  is 
shown  in  Figure  6,  The  waveguide  for  the  structure  used  to  calculate  Figure 
6  is  symmetric  with  respect  to  the  geometry  and  the  refractive  indices.  The 
optical  confinement  factor  is  defined  as  the  overlap  of  the  optical  mode  shape 
with  the  active  region.  That  is, 


f0davUy)dy 

I-ooEo(y)dy 


Eo  (y)dy 


(36) 


13 


3.8 


3.6 

X 

CD 

"O 

c 

CD 

3.4  .> 

-»— » 
o 

CC 


a; 

CL 


3.2 


3 


Figure  6:  Normalized  optical  mode  distribution  and  refractive  index  variation 
with  position  for  a  symmetric  waveguide. 

4.2  Asymmetric  Waveguiding:  SHEDS  Lasers 

Long  wavelength  (Ao  =  1.55  /irn)  laser  structures  utilizing  asymmetric  waveg¬ 
uides  have  been  shown  to  be  advantageous  experimentally  [19j ,  [15j ,  [13J  and 
theoretically  [26j,  [21J ,  |23j,  |22j,  [24J.  Such  lasers  are  known  as  Super  High 
Efficiency  Diode  Structures,  or  SHEDS  lasers.  The  large  density  of  free  carri¬ 
ers  in  the  active  region  of  a  laser  diode  gives  rise  to  free  carrier  absorption  |21j 
and  bimolecular  recombination  [25J  in  the  p-  and  n-doped  core  and  cladding. 
These  phenomena  reduce  the  external  efficiency  of  the  laser.  Free  carrier 
absorption  occurs  as  a  result  of  the  overlap  of  the  optical  mode  shape  with 
the  highly  doped  cladding  layers  [26J .  Quantum  mechanical  calculations  have 
proved  that  the  contributions  to  free  carrier  absorption  by  interconduction 
and  intraconduction  band  absorption  are  negligible  in  laser  diodes  [12J.  The 
dominant  process  is  intervalence  band  absorption  (IVBA)  [22J. 

Proposed  solutions  to  decrease  the  free  carrier  absorption  include  broad¬ 
ening  the  thickness  of  the  core  for  better  optical  confinement,  or  using  an 
asymmetric  waveguide  with  respect  to  the  position  of  the  quantum  well. 
The  effects  of  the  latter  are  explored  in  this  section.  The  p-cladding  is  highly 
lossy,  due  to  two  effects:  the  p-material  has  a  larger  free-carrier  absorption 
cross-section  than  that  of  the  n-material,  and  a  higher  doping  concentration 
to  compensate  for  the  low  hole  mobility  (relative  to  electron  mobility)  [23J. 
Therefore,  optimal  designs  seek  to  reduce  the  overlap  of  the  optical  mode 
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Figure  7:  Normalized  optical  mode  distribution  and  refractive  indices  for  an 
asymmetric  waveguide  (geometry  and  material  properties  from  [22J ) .  Note 
the  overlap  of  the  optical  mode  with  the  n-cladding  and  reduced  overlap  with 
the  p-cladding  relative  to  the  symmetric  waveguide  of  Fig.  6, 


shape  with  this  layer.  This  is  achieved  not  only  by  geometric  asymmetry, 
but  by  index-guiding  the  optical  mode  by  using  a  material  with  higher  refrac¬ 
tive  index  for  the  n-cladding  relative  to  that  of  the  p-cladding.  Index-guiding 
in  this  manner  does  increase  the  overlap  of  the  optical  mode  shape  with  the 
n-cladding.  This  is  illustrated  in  Fig.  7,  which  uses  material  properties  and 
dimensions  from  |22j.  However,  the  n-cladding  can  be  made  much  less  lossy 
than  the  p-cladding  in  terms  of  free  carrier  absorption,  because  the  higher 
electron  mobility  allows  for  a  lower  magnitude  of  donor  doping  concentra¬ 
tion  to  achieve  the  same  electrical  conduction  properties.  The  optical  mode 
shape  determines  the  free  carrier  absorption  at  a  point,  and  therefore  the 
volumetric  heating  term  q'”pt. 

5  Eigenenergy  Calculation 

In  the  quantum  well,  the  conduction  and  valence  bands  develop  into  discrete 
energy  subbands.  The  calculation  of  the  allowed  subband  minima  follows 
a  standard  quantum  mechanical  treatment  (see,  for  example,  [16J,  |3j)  for 
a  single  particle  in  a  finite  potential  well.  The  carriers  are  bound  to  the 
potential  well  by  the  band  energy  offsets  between  the  quantum  well  material 
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and  the  barrier  materials.  The  single  electron,  time-independent  Schrodinger 
equation  is, 


(-|^V2  +  \/)^n  =  «n,  (37) 

where  m  is  the  appropriate  effective  mass  (rri*  for  the  conduction  band,  m*hh 
for  the  heavy  hole  valence  band,  or  m*lh  for  the  light  hole  valence  band),  V 
is  the  form  of  the  potential,  and  if} n  is  the  wave  function  of  the  nth  level. 
The  nth  eigenenergy  is  given  by  En  (for  a  specihc  band  calculation,  En  — > 
Ecn,  Ehhni  or  Eihn ■  We  solve  the  time-independent  Schrodinger  equation  in 
the  dimension  perpendicular  to  the  plane  of  the  quantum  well,  which  we 
have  designated  as  the  ^/-direction.  Assuming  that  the  potential  well  is  of 
the  form, 


V(y) 


0  inside  the  well 
A E  outside  the  well 


(38) 


where  A E  is  the  appropriate  band  offset  (A Ec,  A Ehh,  or  AELh)  between  the 
quantum  well  and  the  barrier  material.  This  parameter  is  calculated  assum¬ 
ing  a  conduction  to  valence  band  offset  ratio  of  3:2  for  AlTGai_,j;As  (barriers) 
to  GaAs  (quantum  well)  hetero junctions  [14J.  The  wave  function  and  its  first 
derivative  must  be  continuous  across  the  heterojunction.  Additionally,  the 
wave  function  must  go  to  zero  far  away  from  the  well.  The  resulting  solution, 
using  (38)  as  the  potential,  is 


tan 


da  yj2m(AE  -  En) 
~2  Ti 


(A  E  \"1/2 

Ur1) 


(39) 


cot 


da  y/2m{AE  -  En) 

2  ti 


/  AE  \ 

- lj  n  =  2,  4,  6,...,  (40) 


where  da  is  the  thickness  of  the  well.  These  transcendental  equations  must 
be  solved  numerically  for  En.  The  solutions  can  represented  graphically  by 
defining  new  parameters, 


V  = 


da  y/2rn  (A E  -  En) 
2  k 


(41) 

(42) 
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Figure  8:  Graphical  representation  of  the  solution  to  Schrodinger’s  equation. 
The  energy  eigenvalues  are  the  intersection  of  the  tan  r)  and  —  cot  ij  curves 
with  the  (  curve.  The  solutions  are  indicated  (in  terms  of  rf)  as  +  marks. 


The  graphical  representation  of  the  solution  is  shown  in  Fig.  8,  The  elec¬ 
tron  eigenenergy  solutions  are  also  plotted  along  with  the  conduction  band 
through  the  well  and  barriers  in  Fig.  9.  The  optical  wavelength  is  inversely 
proportional  to  the  total  difference  in  energy  between  the  lowest  conduction 
and  valence  sub-bands.  Therefore,  this  calculation  must  be  performed  prior 
to  the  optical  mode  and  gain  calculations. 

ft  is  possible  to  couple  the  Poisson,  Schroclinger,  and  current  continuity 
equations  |8j,  [1  lj  in  a  semi-classical  treatment,  but  this  would  require  solving 
Schrodinger’s  equation  at  each  iteration  of  the  drift-diffusion  solver.  There¬ 
fore,  following  [14J,  the  eigenenergies  are  calculated  prior  to  the  drift-diffusion 
module  and  are  assumed  to  be  decoupled. 

6  Conclusion 

A  two-dimensional  electrothermal  model  for  high-powered  diode  lasers  has 
been  developed  and  discussed.  A  drift- diffusion  model  has  been  formulated 
to  include  recombination  mechanisms  relevant  to  a  diode  laser.  A  thermal 
model  coupled  to  the  drift-diffusion  solver  has  been  implemented  to  model 
the  effects  of  temperature  on  electron-hole  transport.  An  optical  model  has 
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Figure  9:  The  allowed  electron  eigenenergies  in  an  8  nm  GaAs  quantum  well 
(red)  along  with  the  conduction  band  profile  through  the  y-direction. 


been  developed  to  determine  the  optical  confinement  factor,  and  therefore 
the  portion  of  light  absorbed  by  the  waveguiding  layers.  Initial  results  for 
an  asymmetric  waveguide  system  have  shown  that  the  optical  mode  can  be 
shifted  into  the  n-doped  cladding  layer,  which  is  less  lossy  than  the  p-doped 
cladding.  Further,  the  coupling  of  the  optical  gain  to  the  complex  index  of 
refraction  used  in  the  optical  mode  calculation  has  been  discussed,  as  well  as 
the  relation  of  optical  gain  to  the  stimulated  emission  rate  used  in  the  drift- 
diffusion  model.  Lastly,  a  decoupled,  one-dimensional  quantum  mechanical 
simulation  has  been  solved  to  determine  the  energy  states  in  the  quantum 
well. 
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Nomenclature 


Mqw 


Oil 


AEf 

A 

r 

h 

K 

'Vj 

C 

[Jjn 

flp 

V-.  V 


PredJ 

Tn 

Tp 

7~in 

e0 

in 

ip 

esi  e 


Quantum  well  transition  matrix  element 
Internal  losses 

Quasi-Fermi  level  separation 
Spin-orbit  splitting  energy 
Optical  confinement  factor 
Reduced  Planck  constant,  ft  =  h/( 2n) 

Thermal  conductivity 
Wavelength  in  vacuum 
Lorentzian  broadening  function 
Electron  mobility 
Hole  mobility 

Divergence  operator,  gradient  operator,  respectively 

Electrostatic  potential 

Equilibrium  potential  distribution 

Reduced  density  of  states  for  jth  conduction  band 

Shockley-Read-Hall  relaxation  time  for  electrons 

Shockley-Read-Hall  relaxation  time  for  holes 

Intraband  relaxation  time 

Optical  mode  shape 

Electron  current  density 

Hole  current  density 

Relative  static  permittivity,  vacuum  permittivity,  and  abso¬ 
lute  static  permittivity,  respectively 

Electron  quasi-Fermi  potential 
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<fip 

Co 

Ec 

Ef 

Eg 

Ei 

Ev 

ECn 

EFc 

EFv 

Ejn 

fcn,j 

fvn,j 

9 

J 

k0 

kF 

L 

m* 

m* 

ml 


M, 


avg 


m 


hh 


m 


lh 


m. 


red,  j 


Hole  quasi-Fermi  potential 
Speed  of  light  in  vacuum 
Conduction  band  minimum  (bulk) 

Equilibrium  Fermi  energy 
QW  band  gap 
Photon  energy 

Valence  band  minimum  (bulk) 

nth  quantized  conduction  band  minimum 

Electron  quasi- Fermi  energy 

Hole  quasi-Fermi  energy 

nth  quantized  j  valence  band  minimum 

Modified  Fermi  distribution  for  conduction  band 

Modified  Fermi  distribution  for  valence  band 

Optical  gain 

j  =  hh,  lh  for  heavy  holes  and  light  holes,  respectively 
Wave  number  in  vacuum 
Boltzmann  constant 

Cavity  length,  perpendicular  to  the  mirror  planes 
Modified  effective  mass 

Conduction  band  density  of  states  effective  mass 
Valence  band  density  of  states  effective  mass 
Bulk  transition  matrix  element 
Heavy  hole  effective  mass 
Light  hole  effective  mass 

Reduced  effective  mass  for  the  jth  valence  band 
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n 


Electron  concentration 


N,  N+,  NJ 
Nc 

rii 

nr,  nr 
Nv 


Tl g,eS 


nm 

Nph 


P 

Ppo 

q 

Cl' Joule 

R 

R\i  R-2 

Kp 

Rst 

T 

Ths 

V 

va 


Doping  profile  as  a  function  of  position,  ionized  donor  con¬ 
centration,  ionized  acceptor  concentration,  respectively 

Effective  density  of  states  for  conduction  band 

Intrinsic  carrier  density 

Complex  refractive  index  and  its  real  part,  respectively 
Effective  density  of  states  for  valence  band 
Group  effective  refractive  index 

Equilibrium  electron  concentration  in  the  n-doped  region 
Photon  density 

Complex  effective  refractive  index  eigenvalue 
Population  inversion  factor 
Hole  concentration 

Equilibrium  hole  concentration  in  the  p-doped  region 

Electron  charge 

Volumetric  Joule  heating 

Recombination  mechanisms 

Mirror  reflectivites 

Spontaneous  emission  rate  emitted  into  the  lasing  mode 

Stimulated  emission  recombination  rate 

Temperature 

Heat  sink  temperature 

Volume  of  the  active  region 

Applied  voltage 

Group  velocity  of  photons 
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